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ABSTRACT 


Inventory  systems  with  returns  are  systems  in  which  there  are 
units  returned  in  a  repairable  state,  as  well  as  demands  for  units  in  a 
serviceable  state,  where  the  return  and  demand  processes  are  independent. 
We  consider  the  inventory  control  of  a  single  item  at  a  single  location 
in  which  the  stationary  return  rate  is  less  than  the  stationary  demand 
rate.  This  necessitates  an  additional  occasional  procurement  of  units 
from  an  outside  source.  The  objectives  of  this  paper  are  to  develop  a 
cost  model  of  this  system  managed  under  a  continuous  review  procurement 
policy,  and  to  develop  a  solution  method  for  finding  the  policy 
parameter  values.  The  key  to  the  analysis  is  the  use  of  a  normally 
distributed  random  variable  to  approximate  the  steady-state  distribution 
of  net  inventory. 


1.  Introduction 


Many  models  have  been  developed  during  the  past  15  years  pertaining 
to  various  aspects  of  managing  repairable  item  inventory  systems  (e.g., 
[l]»[h],[10],[ll],[12],[15]  and  [16]).^  Most  of  these  models  contain  the 
assumption  that  the  failure  of  a  unit  simultaneously  generates  a  demand 
for  a  unit  of  exactly  the  same  type,  i.e.  the  demand  process  for  service¬ 
able  units  and  the  return  processes  of  failed  units  are  perfectly  correla¬ 
ted. 

In  certain  instances,  however,  this  assumption  of  perfect  correlation 
between  the  demand  and  return  processes  is  not  valid.  For  example,  this 
can  occur  in  situations  where  equipment  is  leased,  rented,  and/or  sold, 
such  as  found  in  the  telephone,  computer  and  copying  machine  industries. 
Returns  do  not  necessarily  correspond  to  failures  in  these  cases,  but 
rather  to  lease  or  rental  expirations.  At  the  time  a  unit  is  returned,  it 
may  have  to  go  through  a  repair  or  overhaul  process  before  reissue.  There 
is  no  reason  to  assume  that  with  the  expiration  of  a  lease  or  rental 
agreement  that  the  customer  will  request  a  unit  of  exactly  the  same  type. 
Similarly,  when  a  customer  requests  a  particular  type  of  unit,  there  is 
no  reason  to  assume  that  the  customer  will  return  one  of  exactly  the  same 
type. 

Although  the  assumption  of  partial  dependence  between  the  return 
and  demand  processes  may,  in  general,  be  the  most  realistic,  we  will 
assume  in  our  analysis  that  these  processes  are  independent.  This 
assumption  was  tested  and  found  to  be  valid  for  the  inventory  system  the 


A  repairable  item  is  an  item  which  fails,  but  which  can  be  repaired  and 
subsequently  made  available  to  satisfy  a  future  demand  or  an  existing 
backorder. 
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authors  examined.  We  will  call  inventory  systems  in  which  the  return 
and  demand  processes  for  repairable  items  are  independent,  inventory 
systems  with  returns.  It  is  this  type  of  inventory  system  that  will  be 
analyzed  in  detail  in  this  paper. 

This  section  continues  with  a  description  of  the  inventory  system 
we  will  study  and  concludes  with  a  brief  survey  of  the  existing  litera¬ 
ture  on  this  problem.  In  Section  2  we  develop  the  stationary  distribution 
of  two  key  random  variables  that  describe  the  probabilistic  behavior  of 
the  inventory  system.  This  analysis  is  used  as  the  basis  for  a  cost 
model  presented  in  Section  3.  In  Section  4,  we  examine  an  important 
special  case  of  our  inventory  system,  and  in  Section  5,  we  conclude  with 
a  brief  summary  and  some  final  comments. 

The  system  we  will  study  consists  of  a  single  type  of  item  managed 
at  a  single  location.  A  schematic  representation  of  this  inventory  system 
is  given  by  Figure  1.  As  shown,  this  location  is  assumed  to  contain  both 
a  repair  facility  for  returned  units,  and  a  warehouse,  or  storage  facility, 
for  serviceable  inventory. 


Figure  1.  A  Schematic  Representation  of  the  Inventory  System 
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We  assume  returns  of  repairable  units  occur  as  a  Poisson  process  with 
rate  Y,  and  demands  for  serviceable  units  occur  as  a  Poisson  process  with 
rate  X.  These  two  processes  are  also  assumed  to  be  independent.  Y  is 
assumed  to  be  less  than  X,  so  that  an  occasional  procurement  of  units 
from  an  outside  source  is  required.  Units  procured  in  this  manner  arrive 
in  a  serviceable  state  x  time  units  after  they  are  ordered. 

The  repair  facility  behaves  as  a  first-come,  first-served  queueing  sys¬ 
tem  with  Poisson  arrivals  (the  Poisson  returns).  All  returned  units  require 
repair,  and  repair  times  of  returned  units  are  independent.  Since  Y  <  X, 
the  repair  system  is  always  operating  as  long  as  repairables  are  present. 

No  other  assumptions  about  the  queueing  repair  system  (e.g.  service  time  dis¬ 
tribution,  number  of  repair  servers)  are  made. 

The  output  of  this  queueing  repair  system  is  input  to  the  stock  of  on- 
hand  serviceable  inventory,  as  is  the  arrival  of  outside  procurement  orders. 

All  demands  that  are  not  immediately  satisfied  are  assumed  to  be  back¬ 
ordered  . 

We  define  'net  inventory'  at  a  point  in  time  to  be  the  number  of  on-hand 
serviceable  units  in  the  storage  facility  minus  the  number  of  outstanding 
backorders.  We  also  define  'inventory  position'  at  a  point  in  time  to  be  the 
sum  of  net  inventory,  the  number  of  units  in  the  repair  queueing  system,  and 
the  number  of  units  on-order  from  the  outside  procurement  source. 

Let  I(t)  =  the  inventory  position  at  time  t, 

N(t)  =  the  net  inventory  at  time  t, 

R(t)  =  the  number  of  units  in  the  repair  queueing  system  at  time  t, 
P(t)  =  the  number  of  units  on-order  from  the  outside  supplier  at 


time  t 
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0(t)  =  the  on-hand  serviceable  inventory  at  time  t, 
and  B(t)  =  the  number  of  outstanding  backorders  at  time  t. 

Then 

I(t)  =  N(t)  +  R(t )  +  P(t), 
and  N(t)  =  0(t)  -  B(t). 

Our  final  assumption  concerns  the  form  of  the  procurement  policy. 

We  assume  that  a  continuous  review  (Q,r)  procurement  policy  is  followed, 
i.e.  when  inventory  position  drops  below  r  +  1,  an  order  for  Q  units 
is  immediately  placed.  Since  the  repair  queueing  system  is  assumed  to  be 
continuously  operating,  our  objective  is  simply  to  find  values  of  Q  and 
r. 

The  objectives  of  this  paper  are  to  develop  a  model  of  this  single 
item,  single  location  inventory  system  with  returns,  and  to  present  a 
procedure  for  finding  the  values  of  the  policy  parameters. 

Only  a  few  papers  have  been  published  on  inventory  systems  with  returns. 
These  papers  contain  simplifying  assumptions  which  make  them  of  limited 
practical  value.  Heyman  [6,7]  considers  optimal  disposal  policies  for  a 
single-item  inventory  system  with  returns;  but  his  assumptions  include 
instantaneous  outside  procurement  (implying  no  backorders  or  lost  sales) 
and  no  fixed  cost  of  ordering  (implying  no  lot-size  reordering).  Hoadley 
and  Heyman  [8]  consider  a  two-echelon  inventory  system  with  outside  pro- 
curment,  returns,  disposals,  and  transshipment;  but  their  model  is  a  one- 
period  model,  and  all  of  the  mentioned  transactions  are  assumed  to  occur 
instantaneously.  Simpson  [16]  develops  the  optimum  solution  structure 
for  a  finite-horizon,  periodic  review  model  of  an  inventory  system  with 
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returns.  His  model  allows  for  correlation  between  the  return  and  demand 
processes.  Backlogging  is  permitted,  but  both  repairs  and  outside  pro¬ 
curements  are  assumed  to  be  instantaneous. 

For  the  most  part,  the  methods  of  analysis  in  these  three  papers  rely 
heavily  upon  the  assumptions  of  instantaneous  repair  and  procurement. 

Their  approaches  are  of  little  use  when  analyzing  situations  in  which  re¬ 
pair  and  procurement  times  are  not  zero. 

Finally,  Schrady  [14]  solves  for  repair  carcass  and  procurement  lot- 
sizes  for  a  completely  deterministic  system.  Gajdalo  [2]  extends  this  to 
a  'continuous  review  repair  policy'  for  an  inventory  system  with  stochastic 
(compound  Poisson)  returns  and  demands.  He  uses  computer  simulation  to 
test  several  heuristics  for  computing  the  reorder  point  and  reorder  lot- 
sizes  for  both  procuring  and  repairing  items.  All  lead-times,  including 
repair  times,  are  assumed  constant. 

2.  Analysis 

The  analysis  begins  with  the  derivation  of  the  steady-state  distribu¬ 
tion  of  inventory  position.  This  result  is  used  in  the  derivation  of  an 
approximation  of  the  steady-state  distribution  of  net  inventory,  and  is 
followed  by  a  discussion  of  the  accuracy  of  the  approximation. 

2 . 1  Derivation  of  the  Stationary  Distribution  of  Inventory  Position 

Changes  in  the  state  of  the  inventory  position  are  caused  only  by 
demands  and  returns.  State  i  (i  =  r+l,r+2,...)  can  be  entered  from 
state  i  +  1  when  a  demand  for  a  serviceable  item  occurs;  state 
j  (j  =  r+2,r+3,...)  can  be  entered  from  state  j  -  1  when  an  item  is 
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returned.  In  addition,  state  r  +  y  can  also  be  reached  from  state 
r  +  1  when  a  serviceable  item  is  demanded  (an  order  for  Q  units  is 
placed  immediately  when  the  inventory  position  drops  below  r  +  1). 

The  time  between  state  transitions  is  exponentially  distributed,  since 
the  return  and  demand  processes  are  Poisson  processes.  The  state  tran¬ 
sition  flow  diagram  is  given  in  Figure  2,  with  the  transition  rates  as 
indicated. 

A 


Figure  2.  State  Transition  Flow  Diagram  For  Inventory  Position 

Let  u.  =  lim  Prob(I(t)  =  r  +  1  +  i),  the  stationary  probability 
that  inventory  position  is  equal  to  r  +  1  +  i.  This  limit  exists  because 
the  states  of  this  system  are  the  states  of  an  irreducible,  ergodic, 

Markov  chain  [13].  The  steady-state  balance  equations  correponding  to 
this  system  are 


(A+y)uq  =  +  Aux, 

(A+Y)u^  =  Yl^  +  Au2, 

(A+Y)u2  =  YUl  +  A^, 


(A+Y)uq_i  =  Yuq_2  t  Auq  t  Auo, 
(A+Y)uQ  =  +  AuQ+1, 


(1) 
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and  the  mean  and  variance  of  the  stationary  distribution  of  inventory  posi¬ 
tion  are  given  by 

(4)  E[lim  I(t)]  =  r  +  1  +  G’(l)  =  r  +  l+  ^-  +  ^-, 

fKo  “  ' 

and 

2 

(5)  Var[lim  I(t)]  =  G"(l)  +  G'(l)  -  [G'(l)]2  =  +  ■  ■  —  , 

t-»»  12  (X-Y) 

respectively. 

If  Q  =  1,  Figure  2  is  the  transition  flow  diagram  for  an  M/M/1  queue¬ 
ing  system  in  which  the  ’arrival'  rate  is  y,  and  the  'service'  rate  is 
X.  In  this  case  (3)  reduces  to  the  geometric  distribution,  which  is  the 
well-known  steady  state  distribution  of  the  number  of  customers  present  in 
an  M/M/1  system. 


Note  that  when  y  =  0,  (4),  (5),  and  (3)  reduce  to  the  mean  variance, 
and  probability  distribution,  respectively,  of  a  uniformly  distributed 
random  variable,  a  well  known  result  (see  reference  5). 


2 . 2  An  Approximation  to  the  Stationary  Distribution  for  Net  Inventory 

Next,  we  develop  an  approximation  to  the  stationary  distribution  of 
net  inventory,  which  is  the  basis  for  the  cost  model  used  to  determine 
optimal  values  of  Q  and  r. 

Recall  that  t,  the  procurement  lead-time,  is  constant.  Thus,  any 


units  on-order  at  time  t  -  x  will  have  arrived  by  time  t.  Similarly, 


any  order  placed  after  time  t  -  x  will  not  have  arrived  by  time  t. 


Therefore,  we  see  that 


(6)  N(t )  =  I(t-x)  -  R(t-x)  +  Z(t-x,  t)  -  D(t-x,  t), 


where  R(t-x)  =  the  number  of  units  in  the  repair  system  at  time 
t  "  T, 

Z ( t — x ,  t)  =  the  output  of  the  repair  system  in  the  interval 
(t-T,  t], 

and  D(t-x,  t)  =  the  number  of  demands  in  the  interval  (t-x,  t]. 


R(t-x)  is  subtracted  from  I(t-x)  so  that  we  do  not  double  count  the 
units  in  the  repair  system  at  time  t-x  that  complete  service  by  time 
t.  Therefore,  net  inventory  at  time  t  consists  of  units  on-order  or 
already  serviceable  at  time  t-x  (both  measured  in  I(t-x)),  plus 
those  units  completing  repair  by  time  t-x,  minus  demands  over  the 
interval  (t-x,  t]. 
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Let  us  separately  examine  the  individual  terms  of  (6).  The  steady- 
state  distribution  of  I(t-x)  has  already  been  obtained.  The  number  of 
demands  over  the  interval  (t-T,  t]  is  Poisson  distributed  with  mean 
YT  and  is  independent  of  the  other  three  random  variables  on  the  right- 
hand  side  of  equation  (6). 

The  distributions  of  R(t-x)  and  Z(t-x,  t)  are  readily  available 
for  many  queueing  systems;  but,  they  are  not  independent  of  each  other 
or  of  I(t-t).  The  number  in  the  repair  system  at  time  t  -  t,  R(t-t), 
clearly  depends  on  the  inventory  position  at  time  t  -  t,  I(t-r).  The 
output  of  the  queueing  system  in  (t-t,  t  ],  Z(t-x,  t),  depends  on  the 
number  in  the  repair  system  at  the  start  of  the  interval,  R(t-x).  This 
latter  dependence  decreases  as  x  increases,  and  one  could  assume  as  an 
approximation  that  R(t-x)  and  Z(t-x,  t)  are  independent  for  values  of 
T  that  are  large  relative  to  the  mean  repair  time.  However,  the  dependence 
between  R(t-x)  and  l(t-x)  cannot  be  ignored.  The  joint  distribution 
of  R(t-x)  and  I(t-x)  is  difficult  to  develop  analytically.  Consequently, 
an  approximation  to  the  distribution  of  net  inventory  will  be  developed, 
using  (6),  rather  than  developing  the  exact  distributions. 

The  normal  distribution  was  chosen  to  be  a  continuous  approximation 
to  the  steady-state  distribution  of  net  inventory.  This  was  done  for  two 
reasons.  First,  we  observed  that  the  steady  state  distribution  of  net 
inventory  for  several  test  cases  (obtained  via  simulation)  resembled  a 
normal  distribution.  Second,  a  theoretical  basis  exists  for  the  use  of  the 
normal  distribution.  Note  that  when  the  effect  of  the  return  and  repair 
processes  are  small,  I (t-t)  is  close  to  being  uniformly  distributed.  Also, 
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the  Poisson  random  variable  D(t-T,  t)  is  well  approximated  by  a  normally 
distributed  random  variable  if  At  >  10.  Then  N(t)  is  close  to  being 
the  difference  between  a  uniform  random  variable  and  a  normally  distributed 
random  variable,  which  itself  is  close  to  being  a  normally  distributed 
random  variable . 

Equation  (6)  is  used  to  determine  the  mean  u  and  to  approximate  the 
2 

variance,  o  ,  of  this  normal  distribution.  Letting  t  -*■  00 ,  we  have 

V  =  E(N(t) )  =  E(I(t-T))  -  E(R( t- t) ) 

(7)  +  E( Z( t-T,  t))  -  E(D(t-T,  t)) 

=  r  +  1  +  pp  +  p—  -  E(R(t-T ) )  +  yt  -  At  , 

using  (5),  and  noting  that  the  expected  output  of  a  queueing  system  over 

an  interval  is  equal  to  the  expected  input  over  an  interval  of  the  same 

2 

length.  Also,  by  ignoring  covariance  terms  we  approximate  o  by 

o2  =  Var(N(t))  %  Var(.T(t-T))  +  Var(R(t-T)) 

(8)  +  Var(Z(t-T,  t))  +  Var(D(t-T,  t)) 

2 

=  p-p-  +  — 0  +  Var(R(t-T))  +  Var(Z(t-T,  t))  +  At, 

(A-y) 

using  (5).  Note  that  exact  expressions  and  good  approximations  for 
E(R(t-T)),  Var(R(t-T)),  and  Var(Z(t-T,  t))  are  available  for  many 
queueing  systems  (e.g.  see  [3]). 

The  accuracy  of  the  normal  approximation,  whose  mean  and  variance 
are  given  by  (7)  and  (8),  was  tested  using  an  incomplete  factorial  exper¬ 
iment.  The  variable  factors  were  the  number  of  repair  servers,  the 
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repair  service  distribution,  the  repair  system  traffic  intensity,  the 
procurement  lead-time,  t,  the  procurement  lot-size  Q,  and  the  ratio 
Y/X.  In  each  test  case,  the  accuracy  of  the  normal  approximation,  was 
measured  by  finding  the  area  between  the  normal  curve  and  the  curve  re¬ 
presenting  the  continuous  version  of  the  distrioution  of  net  inventory 
which  was  obtained  via  simulation. 

The  conclusion  drawn  from  this  experiment  was  that  the  major  factor 
affecting  the  accuracy  of  the  normal  approximation  is  the  ratio  of  the 
return  rate  to  the  demand  rate,  Y/X.  In  fact,  the  normal  approximation 
is  quite  accurate  when  Y/A  <  .6.  However,  a  visual  inspection  of  the 
normal  curves  revealed  that  the  normal  approximation  was  a  good  one  in  the 
left-hand  tail  of  the  distribution  of  net  inventory  in  all  the  test  cases. 
As  in  most  inventory  problems,  the  tail  of  the  distribution  is  all  that 
is  needed  to  determine  optimal  values  for  Q  and  r.  (The  reason  for 
this  will  be  discussed  in  Section  3.)  Thus  the  normal  approximation  was 
deemed  to  be  an  acceptable  approximation  for  all  realistic  situations. 

3.  Cost  Model  and  Solution  Method 

The  optimization  model  we  wish  to  study  includes  a  fixed  procurement 
order  cost,  a  holding  cost,  and  a  time-weighted  backorder  cost.  In  partic¬ 
ular,  let 


A  =  the  fixed  procurement  order  cost  ($ /procurement  order), 
h  =  the  holding  cost  ($/unit-year) , 
and  it  =  the  backorder  cost  ($/unit-year) . 

Our  objective  function  K,  will  be  the  expected  sum  of  annual  procurement 
ordering  costs , holding  costs,  and  backorder  costs.  It  will  be  evaluated  by 
taking  the  sum  of 


1)  Ax  (the  expected  number  of  orders  placed  per  year), 

2)  h  x  (the  expected  serviceable  on-hand  inventory  at  a  random  point 


in  time) , 

and  3)  it  x  (the  expected  number  of  outstanding  backorders  at  a  random 
point  in  time). 

Both  the  expected  on-hand  inventory  and  expected  backorders  at  a  random 
point  in  time  will  be  calculated  using  the  normal  approximation  to  the  dis¬ 
tribution  of  net  inventory. 

Note  that  we  need  not  consider  holding  costs  charged  against  units 
in  repair.  Due  to  the  assumption  that  no  inserted  idleness  in  the 
queueing  repair  system  is  allowed,  these  holding  costs  are  independent 
of  the  values  of  the  procurement  policy  parameters. 

Let  <p(  • )  and  4>(  • )  be  the  standard  normal  density  and  standard 
normal  cumulative  functions,  respectively;  i.e.  let 

o 

X 

4>(x)  =  -  •  e  ,  and 

x 

$(x)  =  /  <j>(t)dt. 

t=-  OO 

Let  h(x)  be  the  normal  density,  which  is  the  continuous  approxima¬ 
tion  to  the  steady-state  distribution  of  net  inventory,  and  which  has  mean 
2 

P  and  variance  a  given  by  (7)  and  (8),  respectively.  Thus  the  expected 
number  of  backorders  at  any  point  in  time  is 


(9) 


o<M-)  -  u®(-  -), 

o  a 


as  is  derived  in  reference  5. 

Since 

E(on-hand  inventory)  =  E( inventory  position)  +  E(backorders ) 

-  E(number  in  repair)  -  ECnumber  on-order), 

we  have  the  expected  on-hand  inventory  equal  to 

(10)  r  +  1  +  -2—  +  +  o^(-— )  -  p®(-  -)  -  E(R(t) )  -  (X-y)t. 

i  a-y  o  o 

Note  that  the  last  term,  the  expected  amount  on-order  at  any  point  in  time, 
is  equal  to  the  rate  at  which  demands  are  ultimately  met  by  outside  procure¬ 
ment,  A  -  y,  times  the  constant  procurement  lead-time,  x. 

2 

In  what  follows,  it  will  be  easier  to  think  of  u  and  a  as  functions 
of  r  and  Q.  Specifically,  let 


(11) 

Q 

P  =  r  +  —  +  c 

and 

2 

(12) 

2  Q 

°  =  12  +d’ 

where 

(13) 

c  =  \  -  E(R(t))  -  (A-y)t 

and 

(14) 

d  =  -  -y-0  -  ~  +  Var(R(t) )  +  (A+y)t 

(A-y) 


Finally,  the  rate  at  which  demands  are  met  by  outside  procurement, 
A-Y,  divided  by  Q,  the  procurement  lot-size,  gives  the  expected  number 
of  procurement  orders  placed  per  year. 

Combining  our  previous  results,  we  see  that  the  objective  function 
can  be  expressed  as 
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K  =  U'Q-}-  ♦  -  y<K-  J)] 

+  h[r  +  +  o4>(^-)  -  y4>(-  ^-)  +  c] 

=  L^rllL  +  ( n+h)[o<))(— )  -  y<t>(-  -)]  +  h(r  +  -k?  +  c), 
Q  0  0  2 


where  c  is  given  by  (13). 

This  objective  function  is  not  convex  in  Q,  but  is  convex  in  r. 

This  is  easily  proven  by  showing  that  the  backorder  function  04>(^-)  -  ^-) 

is  convex  in  y.  Since  y=r+-^  +  c,  a  function  that  ie  convex  in 


y  is  also  convex  in  r.  Setting 


0,  we  see  that 


*(-  £)  = 

o  n+h 


Recall,  from  equations  (11)  and  (12),  that  the  mean  y  of  our  ap¬ 
proximating  normal  density  to  net  inventory  is  a  function  of  r,  while 
2 

the  variance,  o  ,  is  not.  Thus,  for  a  fixed  value  of  Q,  the  variance 
of  the  normal  distribution  representing  net  inventory  is  fixed.  Only  the 
mean,  or  'location'  of  the  curve,  is  decided  by  choosing  a  value  of  r. 
Therefore,  equation  (16)  indicates  that  once  the  variance  is  fixed,  the 
'location'  of  the  normal  curve  should  be  chosen  so  that  the  cumulative 
area  to  the  left  of  the  y-axis  is  Jj^-,  as  illustrated  in  Figure  3. 


shaded  area  = 


Figure  3.  Location  of  the  Formal  Curve. 
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In  most  real  situations,  the  backorder  cost  n  is  large  compared  to 
the  holding  cost  h.  This  makes  the  fraction  small.  Recall  that 

this  fraction  is  the  area  to  the  left  of  the  y-axis  under  the  normal  curve. 
The  expected  number  of  backorders  is  calculated  using  equation  (9),  and  the 
expected  on-hand  inventory  is  calculated  in  equation  (10)  also  using  (9). 
Thus,  accuracy  of  the  normal  approximation  is  required  only  in  the  tail  of 
the  distribution,  since  is  usually  small. 

Returning  to  (16)  and  rewriting  it  in  terms  of  r*  and  Q*,  the  opti¬ 
mal  values  of  Q  and  r,  respectively,  we  have 


For  a  fixed  value  of  Q,  the  optimal  value  of  r  is  given  by  equation 
(17). 

To  find  the  optimal  value  of  Q,  one  can  rewrite  equation  (15)  in  terms 
of  r  and  Q.  Using  equation  017')  to  write  the  objective  function  solely 
as  a  function  of  Q,  (15)  simplifies  to 


(18) 
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This  can  be  seen  to  be  a  convex  function  of  Q.  While  the  original 
objective  function,  K,  is  not  convex  everywhere  in  both  Q  and  r, 
upon  deriving  an  optimality  condition  (17),  K  is  convex  in  both  Q 
and  r  over  the  region  of  interest.  Setting  ^  =0,  we  find  that  Q* 
is  the  value  of  Q  that  satisfies 


(19) 


where 


12(A-y)A 


(20) 


a  =  (  Ti-*-h  ) d> (  ^(■"—r-)  ) 

n+h 


A  Fibonacci  search  may  be  used  to  find  0*  since  the  function  on 
the  left  side  of  (19)  increases  with  Q.  (Note,  also,  the  similarity  to 
the  usual  lot-size  formula.  Ignoring  some  of  the  constants,  (19)  is  rough¬ 
ly  of  the  form 

Q  =  •  constant  .  ) 

This  relation  is  independent  of  r.  Thus,  once  Q*  is  found,  r*  is  found 
using  (17). 

We  conclude  this  section  with  a  numerical  example.  Suppose  that  demands 
for  an  item  occur  as  a  Poisson  process  with  rate  A  =  600/year.  The  cost  of 
holding  a  serviceable  unit  in  on-hand  inventory  is  h  =  $200/year;  the  back¬ 
order  cost  is  n  =  $800/backordered  unit  per  year;  the  cost  of  placing  an 
order  is  A  =  $1000/order;  the  procurement  lead-time  is  t  =  .1  years;  returns 
occur  as  a  Poisson  process,  with  rate  y  =  500/year,  to  a  single-server  repair 
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system  with  an  exponential  service  time  distribution  with  rate  n  =  600  units 
per  year. 

The  solution  procedure  starts  by  calculating  the  expected  number  in  re¬ 
pair  and  the  variance  of  the  number  in  repair.  For  an  M/M/1  system,  the  ex¬ 
pected  number  in  repair  is  — =  5  and  the  variance  of  the  number  in 

n- t 

repair  is  — =  30. 

(I--)2 

Using  (13)  and  (14),  c  =  -9.5,  d  =  169.9,  and  a  =  280.3.  Then, 

Q*  satisfies 


12( 100 ) ( 1000  ) 
280.31 


4281.0,  or 


to  the  nearest  tenth,  Q*  =  42.6.  From  this  and  (17), 


^(.8) 


+  9.5 


3.3  . 


If  Q*  and  r*  are  required  to  be  integer,  the  total  cost  K  can 

be  evaluated  for  the  rounded  combinations  of  these  values.  In 
this  case,  the  optimal  combination  is  Q*  =  43  and  r*  =  3.  For  these 
values,  K  =  $7356.  Also,  u  =  15.0  and  o  =  18.0,  so  that  the  ex¬ 
pected  backorders  at  any  point  in  time,  o^(^-)  -  U'K-  ^-),  is  equal  to 
2.03.  The  total  cost  of  $7356  consists  of  $2326  for  ordering,  $1624  for 
backorders,  and  $3406  for  holding  serviceable  stock. 
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4 .  The  Special  Case  y  =  0 

A  special  case  of  the  system  described  in  Section  1  is  one  in  which 
Y  =  0.  Without  returns  of  repairable  items,  the  model  simplifies  to  the 
standard  (Q,r)  model  with  Poisson  demands  and  a  constant  procurement 
lead-time.  A  common  method  of  solution  to  this  problem  is  found  in 
Chapter  4  of  reference  5.  This  method,  which  will  be  referred  to  as  the 
'standard'  method  (SM),  will  be  compared  to  the  solution  method  derived 
previously  for  the  case  y  =  0.  The  latter  will  be  referred  to  as  the  nor¬ 
mal  approximation  method  (NAM). 

The  SM  is  based  on  the  assumption  that  the  probability  that  demand  over 
a  lead-time  exceeds  r  +  Q  is  negligible.  This  approximation  results  in 
the  following  simplified  objective  function  which  is  jointly  convex  in  Q 
and  r : 

(21)  K  =  ~  +  h(j  t-i+r-Ax)  +  i(ir+h)3(r) , 
where 

00 

(22)  B(v)  =  l  (u-v-l)£(u;  At) 

u=v+l 


^ ~~2 P(v-1;  At)  -  AtvP(v;  At)  +  p(  V+1 .  At), 


(23) 


and 


P( a , p )  =  l  p(x;  p ) , 
x=a 


p(x;  p)  =  e  MpX/x! 


(24) 


The  costs  h,  n,  and  A  are  defined  as  before. 


i 


The  optimal  value  for  r,  call  it  r*,  is  the  largest  integer 
satisfying 

(25)  rp(r;  Ax)  -  (r-Ax)P(r;  Ax)  >  , 

and  the  optimal  value  for  Q,  call  it  Q*,  is  the  largest  integer  sat¬ 
isfying 


(26) 


Q(Q-l)  <  [A  +  e(r)]> 

h  A 


or  Q  =  1. 


The  SM  is  an  iterative  algorithm.  The  algorithm  normally  starts 
•  u  „  /5 AA 

with  Q  =  /  — .  A  value  for  r  is  found  using  (25),  and  a  new 

value  for  Q  is  found  using  (26).  This  procedure  continues  until  the 
values  of  Q  and  r  do  not  change. 

The  NAM  may  be  applied  simply  by  setting  Y  =  0.  Since  inventory 
position  is  uniformly  distributed  between  r  t  1  and  r  +  Q  when 

Y  =  0,  the  mean  and  variance  of  the  stationary  distribution  for  inventory 

2 

position  are  r  +  ^  and  >  respectively.  The  net  inventory 

at  time  t,  N(t),  is  given  by 

(27)  N(t)  =  I(t-x)  -  D(t-t,  t). 

Then  the  mean  and  variance  of  net  inventory  can  be  expressed  as 

2 

Q  2  0 

p  =  r  +  —  +  c  and  o  =  +  d,  where 


30 


(28) 

and 

(29) 


Furthermore,  the  variance  in  this  case  is  exact.  With  only  these  small 
changes  in  the  constants  c  and  d,  r*  and  Q*  are  found  as  before 
using  (17)  and  (19). 

The  biggest  advantage  of  using  this  normal  approximation  method  with 
Y  =  0  is  that  an  iterative  procedure  is  not  required  to  find  Q*  and  r*. 

A  binary  or  Fibonacci  search  quickly  determines  the  value  of  Q*,  from 
which  r*  is  uniquely  determined. 

A  factorial  experiment  was  run  to  test  the  NAM  versus  the  SM,  and  to 

determine  which  factors  influence  this  comparison.  Five  different  values 

-  3 

for  each  of  the  factors  it,  AA,  and  At  were  used,  resulting  in  5  =  125 

test  cases.  The  values  used  are  listed  in  Table  1. 


Table  1.  Data  for  Factorial  Experiment 


n  AA  At 


100 

1.000 

5 

200 

2,000 

10 

500 

5,000 

25 

1,000 

10,000 

50 

2,000 

20,000 

100 

The  holding  cost  was  fixed  at  h  =  100. 
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Recall  that  the  SM  is  based  on  an  approximation  to  the  objective 
function,  as  is  the  NAM.  Thus,  in  order  to  test  the  solutions  generated 
by  the  two  methods,  we  use  an  objective  function  which  is  the  exact  mea¬ 
sure  of  the  expected  ordering,  holding,  and  backorder  costs.  This  function 
is 

(30)  K  =  ^  +  h(j  +  j  +  r  -  Xt)  +  7i+h)tn(r)  -B(r+Q)3, 

where  B(v)  is  given  by  (22).  (This  expression  is  derived  by  Hadley  and 
Whitin  [5  ].)  Note  that  the  term  ^<Ti+h)[£(r+Q)]  was  dropped  from  the  ob¬ 
jective  function  in  the  SM  in  order  to  simplify  the  analysis. 

The  results  of  the  experiment  may  be  summarized  as  follows.  The  NAM 
yielded  an  expected  annual  total  cost  (using  K  as  defined  in  (30))  less 
than  or  equal  to  that  obtained  by  the  SM  in  110  of  the  125  cases,  or  88% 
of  the  time.  It  yielded  a  value  for  K  which  was  strictly  less  than  that 
obtained  by  the  SM  in  95  cases,  or  77%  of  the  time.  In  these  95  cases,  the 
average  reduction  in  total  cost  was  approximately  1%  with  the  extreme  case 
having  a  reduction  of  10%.  In  the  15  cases  in  which  thr  normal  approxima¬ 
tion  method  yielded  a  higher  value  of  K,  the  average  increase  was  2.5%, 
with  the  extreme  case  having  an  increase  of  6%. 

The  only  factor  which  had  a  noticeable  effect  on  the  difference  of  the 
results  between  the  two  methods  was  Xr,  the  lead-time  demand.  This  is  ex¬ 
plained  by  the  fact  that  the  normal  distribution  becomes  a  better  continuous 
approximation  to  the  Poisson  distribution  as  the  mean  of  the  Poisson  distri¬ 
bution  increases.  Thus,  as  Xt  increases,  the  normal  distribution  becomes 
a  better  continuous  approximation  to  the  stationary  distribution  of  net  in¬ 


ventory. 
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To  illustrate  these  observations,  in  the  25  cases  in  which  It  =  5, 
the  normal  approximation  method  yields  a  higher  value  of  K  in  eight  of 
the  cases,  or  32%  of  the  time.  The  average  increase  in  K  in  these  eight 
cases  is  2%,  and  the  largest  increase  is  6%.  In  the  25  cases  in  which 
Ax  =  10,  the  normal  approximation  method  yields  a  higher  value  of  K  in 
only  five  of  the  cases,  or  20%  of  the  time.  Here,  in  these  five  cases, 
the  average  increase  is  only  1%  with  the  largest  increase  being  3%.  Yet  of 
the  remaining  75  cases,  in  which  we  have  At  >  25,  the  normal  approxima¬ 
tion  method  yielded  a  higher  value  of  K  in  only  2  of  these  cases. 

Recall,  also,  that  the  SM  is  an  iterative  one.  The  average  number  of 
iterations  necessary  to  solve  these  125  test  cases  using  the  SM  was  3.14. 

The  number  of  iterations  ranged  from  two  to  ten,  and  increased  with  an 
increase  in  At.  In  the  50  cases  in  which  At  >  50,  the  average  number 
of  iterations  was  4.18,  and  rose  to  4.90  in  the  25  cases  in  which  At  =  100 

The  difference  in  computation  time  on  an  IBM  370/168  also  increased 
as  At  increased.  The  25  test  cases  in  which  At  =  5  took  .24  seconds 
of  execution  time  to  solve  via  the  NAM,  and  took  .29  seconds  with  the  SM. 

In  the  25  test  cases  in  which  At  =  50,  however,  the  NAM  took  .23  seconds 
while  the  SM  required  .62  seconds. 

In  conclusion,  the  use  of  the  NAM  should  be  determined  by  the  magni¬ 
tude  of  At.  In  the  25  cases  in  which  At  =  5,  the  two  methods  yielded 
virtually  identical  values  of  K  (20  of  the  25  cases  had  values  of  K  with 
in  1%  of  each  other).  Also,  only  two  iterations  were  needed  using  the  SM 
in  24  of  the  25  cases.  Thus,  for  At  small,  the  two  methods  appear  to  be 
about  equal  in  accuracy  and  computational  complexity. 
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However,  as  soon  as  At  exceeds  10,  both  the  accuracy  of  the  NAM 
improves,  and  the  number  of  iterations  necessary  to  solve  for  Q*  and 
r*  using  the  standard  method  increases.  Thus,  the  NAM  is  recommended 
for  this  cost  model  when  At  >  10  for  both  accuracy  and  computational 
reasons . 

5 .  Summary  and  Concluding  Comments 

We  have  shown  how  to  find  near-optimal  parameter  values  for  a  pro¬ 
curement  policy  for  certain  inventory  systems  with  returns.  The  key  was 
the  use  of  a  normal  approximation  to  the  steady-state  distribution  of  net 
inventory.  This  led  to  the  solution  of  a  cost  model  which  was  convex  in 
the  procurement  parameters.  When  y  =  0,  this  solution  method  provides 
an  accurate  and  efficient  alternative  to  the  currently  accepted  iterative 
algorithm  for  the  standard  Poisson  demand,  constant  lead-time  inventory 
problem  without  returns. 

In  our  model,  we  assumed  the  procurement  policy  to  be  a  stationary 
(Q,r)  policy.  This  policy  is  by  no  means  the  optimal  one.  In  reference 
9,  it  is  shown  that,  for  the  special  cases  of  M/M/1  and  M/G/°°  queueing 
repair  systems,  one  can  lower  total  expected  costs  by  redefining  inven¬ 
tory  position,  and  allowing  variable  reorder  points  as  follows.  Inventory 
position  is  redefined  to  be  net  inventory  plus  the  number  of  units  on-order. 

The  analysis  proceeds  exactly  as  described  in  Section  2  (with  some  of  the 

2 

constants  redefined).  This  results  in  a  reduction  in  o  ,  the  variance  of 

net  inventory,  since  the  variance  of  the  number  of  units  in  repair  is 
.  .  2 

no  longer  included  in  a  .  The  reorder  point,  expressed  in  terms  of 
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inventory  position,  is  then  a  function  of  the  number  of  units  in  repair, 
rather  than  a  constant.  Reductions  in  total  expected  costs  can  be  achieved 
by  using  a  state  dependent  reorder  point  when  the  variance  of  the  number 
of  units  in  repair  is  very  large.  A  10%  reduction  in  total  expected  cost 
was  achieved  using  the  variable  reorder  point  policy  in  an  M/M/1  repair 
system  with  traffic  intensity  p  =  499/500.  This  is  an  extreme  case,  how¬ 
ever.  The  average  annual  cost  of  using  the  stationary  (Q,r)  policy  was 
within  1%  of  the  average  annual  cost  obtained  using  the  non- stationary  one 
in  almost  all  test  cases.  Since  this  is  the  case,  and  since  a  stationary 
(Q,r)  policy  is  easy  to  use,  the  stationary  (Q,r)  policy  is  an  attractive 
policy  to  implement. 
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